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METHOD AND SYSTEM FOR IMPROVING 
THE SPATIAL RESOLUTION FOR STRAIN IMAGING 



FIELD OF THE INVENTION 

The present invention relates generally to medical imaging and more specifically to a 
technique for quantifying ultrasound strain imaging. 

BACKGROUND OF THE INVENTION 

Strain imaging has been used in ultrasound medical imaging to differentiate between 
hard tissue and soft tissue. For example, since tumors are generally harder than normal tissues, 
strain imaging may be used evaluate solid breast masses. Strain images are created by 
comparing echo data obtained before and after a slight compression of the tissue. The results of 
the comparison are displayed as an image on which the hard areas appear dark and soft areas 
appear bright. Other black-and-white and other color scale may be employed for display. 
Typically strain imaging is accomplished by comparing the hard and soft tissue in a region of 
interest on a pixel by pixel basis using conventional ultrasound equipment. 

Typically strain images are made from strain measure point by point, pixel by pixel, 
before and after compression. However, there are several problems with this type of method for 
quantifying the strain. First, there are artifacts that can be generated based on the probe usage 
over several points. Second, the system is time-consuming and tedious. Thirdly, other areas 
can distort the measured strain due to shadows, cavities and other areas of anomalies within the 
area being imaged. 

Accordingly, what is desired is a system to more accurately and quickly quantify 
compression-induced strain the traditional brightness mode (B-mode) images. The present 
invention addresses such a need. 

In recent years, several strain-imaging techniques have been proposed and developed by 



various research groups. In order to increase the spatial resolution for strain imaging, the 
adjacent local blocks, which are typically used to determine their displacements should be close 
enough. But there are technical difficulties for accurately determining the difference of 
displacements over these local blocks and generate the poor signal to noise ratio (SNR) for the 
strain images and therefore quantifying the strain measurements. Additionally, the levels of 
echo signals (which are related to the back-scatter coefficients of insonifying tissue) from 
different types of tissues may be quite different, such as the echo signals from healthy breast 
tissue could be significantly different from the echo signals from breast cancer. These 
differences will cause additional errors for our motion estimation. 

Accordingly, what is needed is a system and method for overcoming these problems. 
The present invention addresses such a need. 

SUMMARY OF THE INVENTION 

A method and system for quantification of strain imaging is disclosed. The method 
and system comprises performing a motion analysis on at least two regions of interest. The 
method and system further includes providing a strain estimate for each of the at least two 
regions of interest and comparing the strain estimates of each of the at least two regions to 
quantify the strain for the at least two regions of interest. 

A system and method in accordance with the present invention provides for strain 
quantification based on conventional B-mode images. Using this technique, the strain of 
regions of interest (ROI) defined by users can be determined and quantitative comparisons can 
be effectively made in real time. The strain quantification can be used to determine tissue's 
properties and can potentially be applied in breast imaging as well as cardiac imaging. 



BRIEF DESCRIPTION OF THE DRAWINGS 



Figure 1 is a simplified block diagram of an ultrasound system according to the present 
invention. 

Figure 2 is a flow chart of a system and method in accordance with the present 
invention. 

Figure 3 shows a plurality of blocks generated for a region of interest. 

Figure 4 illustrates two regions of interest in which the strain has been calculated. 

DETAILED DESCRIPTION 

The present invention relates to medical imaging and more specifically to a technique 
for quantifying tissue strain images. The following description is presented to enable one of 
ordinary skill in the art to make and use the invention and is provided in the context of a patent 
application and its requirements. Various modifications to the preferred embodiment and the 
generic principles and features described herein will be readily apparent to those skilled in the 
art. Thus, the present invention is not intended to be limited to the embodiment shown but is to 
be accorded the widest scope consistent with the principles and features described herein. 

Figure 1 is a simplified block diagram of an ultrasound system according to the present 
invention. The user enters the various conventional scan parameters into an input unit 100, 
which typically includes such devices as a keyboard, knobs and buttons. The input unit is 
connected to a processing system 102, which will typically be an electrically connected and 
cooperating group of processors such as microprocessors and digital signal processors; the 
processing system may, however, also be implemented by a single processor as long as it is fast 
enough to handle the various tasks described below. 

As in known systems, the processing system 102 sets, adjusts, and monitors the 
operating parameters of a conventional transmission control circuit 104, which generates and 
applies electrical control and driving signals to an ultrasonic probe 106, which includes an array 



108 of piezoelectric elements. As is well known in the art, the piezoelectric elements generate 
ultrasonic waves when electrical signals of the proper frequency are applied to them. 

By placing the probe 106 against the body of a patient, these ultrasonic waves enter a 
portion 1 10 of the patient's body. By varying the phasing, amplitude, and timing of the driving 
signals, the ultrasonic waves are focused to form a series of scan lines 112 that typically fan out 
from the probe. Several such scan lines are shown extending into the patient's body in Fig. 3. 
A region of interest, that is, the region that the user wants to have an image of, is shown as an 
interrogation region or volume 114. The manner in which ultrasonic scanning signals are 
controlled, generated, and applied to a patient's body is well understood in the art and is 
therefore not described further. The interrogation volume 114 may be scanned using a series of 
substantially adjacent scan planes (each comprising several scan lines) that extend over a known 
depth. 

Returning to Figure 1, ultrasonic echoes from the waves transmitted into the body return 
to the array 108. As is well understood, the piezoelectric elements in the array thereby convert 
the small mechanical vibrations of the echoes into corresponding electrical signals. 
Amplification and other conventional signal conditioning is then applied to the return signals by 
a reception controller 116. This processing includes, as needed, such known signal conditioning 
as time-gating, gain compensation, and diffraction compensation in order to identify the echo 
signals that correspond to each scan plane of the interrogation volume 114. 

The reception controller 116, all or part of which is normally integrated into the 
processing system 102, converts the ultrasonic, radio-frequency (RF) return signals (typically on 
the order of a few to tens of megahertz) into lower frequency ranges for processing, and may 
also include analog-to-digital conversion circuitry. This is well known in the art of ultrasonic 
imaging. The down-converted power values for the two-dimensional interrogation region are 
stored in a memory 1 1 8 as frame data 120, after conventional beamforming. Each set of frame 



data corresponds to one image frame, that is, to a cross section of the interrogation volume. 
Each frame of the image is represented and stored digitally as an array of acoustic power or 
intensity values for the image elements that make up the frame. 

The interrogation region is normally not in the same shape as what the user wants to see 
displayed, and even when it is, the digital acoustic intensity values formed into beams are 
normally not in a form suitable for driving a conventional gray-tone or color display directly. 
The acoustic intensity values for an image frame are therefore applied to a conventional scan 
converter 122, which converts the digital acoustic values into display intensity or brightness 
values that are suitable for driving a display device 124. The display 124 is typically divided 
into a pattern of picture elements or "pixels" that make up an image that the user can view and 
interpret. Scan conversion and display are well-known features of an ultrasonic imaging system 
and are therefore not described further. 

A system and method in accordance with the present invention provides for strain 
quantification based on conventional B-mode images. Using this technique, the strain of 
specific regions of interest (ROI) defined by users can be determined and quantitative 
comparisons can be effectively made in real time. The strain quantification can be used to 
determine tissue's properties and could be potentially applied in breast imaging as well as 
cardiac imaging. 

To more particularly describe the features of the present invention, refer now to the 

following discussion in conjunction with the accompanying figures. In a preferred embodiment, 
a process is utilized to provide the strain quantifying arrangement. This process could be 
implemented in an algorithm which is part of a computer program. The program can be 
implemented in computer readable medium within the processor 1 02. The computer readable 
medium can be implemented in a disk drive floppy drive, CD-ROM, DVD or the like. 



Figure 2 is a flow chart of a system and method in accordance with the present 
invention. First, at least two regions of interest (ROIs) are selected over one of a sequence of 
conventional B-mode images, via step 202. Several small local blocks, which touch the 
boundary of the ROI, are then generated over each of the regions of interest. Figure 3 shows a 
5 plurality of blocks 150 generated for a region of interest 152. The size of the block 150 can not 

be too small, which could cause more statistical uncertainty. For an example, the size of block 
could be selected as 25x25 pixels. The displacement or motion analysis of each block 150 is 
then determined, via step 204. In a preferred embodiment a conventional block matching 
technique is utilized. Accordingly, the block match technique would match a block in the first 
10 rj B-mode image (uncompressed) to a block in the second B-mode image (compressed) by 

pj computing the correlation of the block in the first image with each block of the same size in a 
y \ defined search region of the second image. In an example, the size of a search region could be 
ffl 32x32 pixels. In order to increase the computational speed and accuracy, the position of MS AD 
f== (the minimum of the Sum of Absolute Difference) of image data is searched instead of 
15 p conventional correlation calculation. After determining all displacements for all local small 
C3 blocks surrounding each region of interest, the strain value for each of the ROI can be 
quantitatively determined, via step 206 by the equation: 



ST = 



xl00% 



20 where ai and bi are the displacement components for two blocks, which cross over the boundary 

of a specific ROI, in the direction of i-th A-line. di is a distance between the two blocks, and ii 
and ii are the index of A-line on B-mode image covering that specific ROI. 

In clinical applications, two ways can be used to increase the accuracy of the strain value 
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measurements. One way is to increase the compression from the probe, and another way is to 
increase the size of the regions of interest. For example, an ROI could be selected, which 
contains MxN pixels (M pixels in axial direction and N pixels in lateral direction). The 
correlation length over lateral direction (which is related to lateral B-mode image's spatial 
resolution) is described as n pixels. In that case, the rninumum strain, which can be detected, is 
given by the equation: 

ST . _ 1Q0 ^ q. 



For an example, if the ROI is a square with 5.0x5.0 mm 2 . The B-mode image is formed 
with a 7.5L40 probe with a depth of 4.0 cm. In that case, the minimum displacement estimation 
10 nj is about § = 0.25 pixel and the pixel size is about 0. 1 mm, and lateral resolution (defined as the 
SJ 6-db width of the point spread function) is about 0.5 mm. The minimum strain, which can then 
be detected, based upon the above-identified equation is 

£ ST- 100V2" x0.25 n/ 

? "%n — / — /o — u.z /o 

Oj (5.0/0.1)xV(5.0/0.1)/(0.5/0.1) 

iSSSSK 

Ifthe size of ROI is about 10.0x10.0 mm 2 , the minimum strain, which can be detected, will be 
15 about 0.078%. 

Accordingly, as above mentioned a strain estimate is performed for each ROI based 
upon the motion analysis, via step 206. Thereafter, the strain estimates are compared to quantify 
which region is hardest, via step 208. 

The results from tests show that a method and system in accordance with the present 
20 invention could provide accurate quantitative information on tissue's strain, which could be 

important for quantifying tissue's properties. Figure 4 illustrates two regions of interest, ROI 
202 and ROI 204 in which the strain has been calculated in accordance with the present 
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invention. As is seen, the strain estimate for the ROI 202 is 2.5% and the strain estimate for 
ROI 204 is 1.5%. Accordingly, as is seen ROI 204 has the harder tissue. In clinical applications 
it is desirable to select at least two ROIs to make an adequate comparison for strain 
quantification. The strain for two specific regions of interest (such as ROI 202 and ROI 204) 
can be determined by using the technique over B-mode images which is provided by a probe 
which provides an image depth of 4.0 cm. In this embodiment, the sizes of ROI are about 
10.0x10 mm 2 . 

Although the present invention has been described in accordance with the embodiments 
shown, one of ordinary skill in the art will readily recognize that there could be variations to the 
embodiments and those variations would be within the spirit and scope of the present invention. 
Accordingly, although the analysis has been accomplished through a block match technique 
one of ordinary skill in the art readily recognizes that other types of motion analysis techniques 
could be utilized and they would be within the spirit and scope of the present invention. 
Accordingly, many modifications may be made by one of ordinary skill in the art without 
departing from the spirit and scope of the appended claims. 



